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Abstract 

A generalized form of Wien's displacement law and the blackbody radiation laws of (a) Rayleigh and 
Jeans, (b) Rayleigh, (c) Wien and Paschen, (d) Thiesen and (e) Planck are derived using principles of 
dimensional analysis. This kind of scaling is expressed in a strictly mathematical manner employing 
dimensional 7r-invariants analysis sometimes called Buckingham's 7r-theorem. It is shown that in the 
case of the classical radiation law of Rayleigh and Jeans only one 7r number occurs that has to be 
considered as a non-dimensional universal constant. This n number may be determined theoretically 
or/and empirically. It is also shown that dimensional 7r-invariants analysis yields a generalized form of 
Wien's displacement law. In this instance two 7r numbers generally occur. Consequently, a universal 
function is established that is indispensable to avoid the so-called Rayleigh- Jeans catastrophe in the 
ultraviolet. Unfortunately, such a universal function cannot be inferred from dimensional arguments. It 
has to be derived theoretically or/and empirically, too. It is shown that such a similarity function can 
be deduced on the basis of heuristic principles, when criteria like the maximum condition regarding the 
generalized form of Wien's displacement law, the requirement of the power law of Stefan and Boltzmann, 
and Ehrenfest's arguments regarding the red and the violet requirements are adopted. 

1 Introduction 

Our contribution is focused on the heuristic derivation of blackbody radiation laws, where principles of 
dimensional analysis are considered. The idea on which dimensional analysis is based is very simple. It is 
inferred from the fact that physical laws do not depend on arbitrarily chosen basic units of measurements. In 
recognizing this simple idea, one may conclude that the functions that express physical laws must possess a 
certain fundamental property, which, from a mathematical point of view, is called the generalized homogene- 
ity or symmetry [3]. This property allows the number of arguments in these functions to be reduced, thereby 
making it simpler to obtain them. As Barenblatt [3 J pointed out, this is the entire content of dimensional 
analysis - there is nothing more to it. 

Often, solutions for physical problems, especially in mechanics and fluid mechanics (e.g., [TJ |5J [21 SI |5J 
[TSJ H3] ) and cloud microphysics jjj], can be found on the basis of similarity hypotheses that comprise 

all problem-relevant dimensional quantities and serve to possess the physical mechanisms of these problems. 
Such similarity hypotheses implicitly describe the functional dependence between these dimensional quanti- 
ties in a mathematical form. This does not mean that this functional dependence can explicitly substantiated 
by formulating a similarity hypothesis only. A similarity hypothesis will become successful if a generalize 



2 



Gerhard Kramm and Fritz Herbert 



homogeneity or a symmetry exist. 

If similarity is hypothesized, its mathematical treatment can further be performed by the procedure 
of dimensional 7r-invariants analysis (sometimes also called the Buckingham's 7r-theorem, for details see 
[TJ [21 [31 131 [TDJ During this mathematical treatment the explicit dependence between the problem- 
relevant dimensional quantities is expressed by the non-dimensional 7r-invariants. 

The first attempt to derive a blackbody radiation law on the basis of dimensional consideration was 
performed by Jeans [T^l [H] using the wavelength A, the absolute temperature T, the velocity of light in 
vacuum c, the charge e and the mass m of an electron, the universal gas constant R, and the dielectric 
constant of the ether K. His attempt to derive Wien's [37] displacement law, however, was strongly criticized 
by Ehrenfest [T2J [T3] . Ehrenfest showed that this kind of dimensional analysis which leads to a similarity 
function containing two tt numbers (see Appendix A) was rather arbitrary. In his reply Jeans |21j rejected 
Ehrcnfest's criticisms. Some months later, Ehrenfest |13j also plucked Jeans' additional arguments to pieces. 
Since their debate was mainly focus on dimensional analysis for which a closed theory was not available 
during that time (a first step on this road was made by Buckingham |10j). there was, if at all, only a minor 
interest on this debate [551. 




In the following we will derive a generalized form of Wien's [47 displacement law and various blackbody 
radiation laws from the perspective of dimensional scaling. In doing so, we partly follow the ideas outlined 
by Glaser (as cited by [33]) and Sommcrfeld [33], i.e., we only postulate that similarity exists. In the chapter 
2, the method used in our dimensional scaling is presented in a strict mathematical manner. The application 
of our method in deriving the blackbody radiation laws is described in chapter 3. Here, it is shown that in 
the case of the classical radiation law of Rayleigh [2pJ [30] and Jeans [2UJ only one tt number occurs that has 
to be considered as a non-dimensional universal constant. This tt number may be determined theoretically 
or/and empirically. In the instance of the generalized form of Wien's |47j displacement law that, in principle, 
contains the radiation laws of Wien [35] and Paschen [35] , Rayleigh [2H] , Thiesen [35] and Planck [33 [3H1 13S] 
as special cases, two tt numbers generally occur. Consequently, a universal function (also called the simi- 
larity function) is established that is indispensable to avoid the so-called Rayleigh- Jeans catastrophe in the 
ultraviolet 14J. Unfortunately, such a universal function cannot explicitly be determined by dimensional 
7r-invariants analysis. It has to derive on the basis of theoretical or/and empirical findings, too. It is shown, 
however, that such a similarity function can be deduced on the basis of heuristic principles, when criteria 
like the maximum condition regarding the generalized form of Wien's |47j displacement law, the requirement 
of the power law of Stefan [35] and Boltzmann [7J, and Ehrcnfest's [T3] red and the violet requirements are 
adopted. Since even Planck's radiation law can be derived in such a manner, some historical notes regarding 
the foundation of the quantum theory are briefly gathered in chapter 4. 

2 Dimensional 7r-invariants analysis 

The theoretical foundation of the procedure is linked to various sources, for instance, Kitaigorodskij [23J, 
Barenblatt [TJ ® |3J g] , Herbert [T7J, Pal Arya [35], Brown 0, Sorbjan [33], and Kramm and Herbert [25] 
which are devoted to characteristic scaling problems in fluid dynamics and turbulence, boundary layer mete- 
orology and other physical disciplines. The description mainly follows the guideline of Kramm and Herbert 

Pi- 
Let adopt that, associated with a certain physical problem, we can select a set of characteristic dimension- 
ality quantities, for instance, n variables, parameters or/and constants, Q%, Q2, ■ ■ ■ , Q K that unambiguously 
and evidently represent the arguments of a mathematical relationship. First this "law" is unspecified; there- 
fore it is formally employed as a general postulate, commonly referred to as the similarity hypothesis of the 
problem, which may read 




F(Q 1 ,Q 2 ,...,Q K ) = 



(2.1) 



In its implicit representation Eq. (2.1) declares k — 1 free or independent arguments as well as a transfor- 
mation of the full series of Q j for j = 1 , . . . , k to a series of p non-dimensional invariants TTi for i — 1 , . . . , p 
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in terms of a factorization by powers. Correspondingly, in that mind each ^-expression is defined by 

- Qt QT ■■■(};: \[Q, i i> , (2.2) 

3=1 

and it is necessarily linked with the condition of non-dimensionality 

dim 7Tj = 1 for i = 1, . . . ,p , (2.3) 

where p < k is customarily valid. 

Next, we will suppose that the 7r-invariants can have interdependencies of arbitrary forms, and it may 
exist a corresponding relation 

0(7ri,7r 2 ,...,7Tp) =0 (2.4) 



which is to be understood as an alternative similarity hypothesis to Eq. (2.1 ). In this function the powers 



are basically unknown numbers, and their determination is the proper problem of the so-called Buckingham 



7r-theorem. If there are more than one 7r-invariant, i.e., p > 1, then we have with Eq. (2.4 1 the explicit 
representation 

7Tj = ^(7Ti,7r 2 , . . . ,7T p ) (2.5) 

in which tp may be interpreted as a universal function within the framework of the similarity hypothesis, 



where, according to the implicit formulation (2.4 1, 7r.j (for any arbitrary i e {1, . . . ,p}) is not an argument of 



that universal function. Note that in the special case of p = 1, we will merely obtain one 7r-invariant, that 



is a non-dimensional universal constant. This special case is expressed by Eq. (2.5 1 in the singular form 



7r = const. (2-6) 

(or if = const.). In view to the determination of the powers xj^, we will extend our treatment to the concise 
set of fundamental dimensions, D n for n — l,...,r, such as length L, time T, mass M, temperature 0, 
considering that any quantity's dimension can be analysed in terms of the independent D n by homogeneous 
power factorization. Let that be expressed as 

r 

dim^ D? D? ...I)';. \\ /•>;! forj = l,..., K , (2.7) 

n=l 

in which the powers g n ,j for n — 1, . . . , r and j — 1, . . . , k are known from the relevant quantities Qj ac- 
cording to the hypothesized similarity condition. Note that r < k is valid, where r is the highest number 
of fundamental dimensions that may occur. In other words: for k quantities Qj including r fundamental 
dimensions D n we obtain p = k — r independent non-dimensional invariants, so-called tt numbers. 



Now a straight- forward development of the analytical framework is attained by introducing Eq. (2. 7 1 



together with the factorization by powers from Eq. (2.2) into the condition of non-dimensionality (2.3). In 



doing so, we obtain this basic law as described in the following detailed representation 

k r 

dim 7Ti = IX(II D n' 3 ) ^ = 1 far * = 1, • • • , J> ■ (2.8) 

3=1 n=l 

Combining the two factorizations ■ and Y[ n m this equation enables to rewrite this set of conditions in the 
fully equivalent form 



dim tt 4 = J] dF 3=1 9n ' 3 X3 " =1 for i = 1, . . . ,p . (2.9) 

n=l 
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For the following conclusion, the latter is more suitable than the former. Indeed, we may immediately infer 
from the factorizing analysis in dependence on the bases D n for n = 1, . . . , r, that the set of condition 



' 

3=1 



Xji = for n = 1,.. . ,r and i = 1, . . , ,p , (2-10) 
has to hold since each £>„-exponential factor must satisfy, owing to its mathematical independence, the 



condition of non-dimensionality (see Eqs. (2.3) and ( 2 . 9 1 ) , i.e., to be equal to unity. In matrix notation, Eq 



(2.101 may be expressed by 




{0} 



(2.11) 



dimensional matrix 
G = {g r , K } 



matrix of powers 
A = 



where the notation {0} is an r x p matrix, and each column of the matrix of powers, A, is forming so-called 
solution vectors Xj for the invariants Hi for i = 1, . . . ,p. The set of equations (2.11 1 serves to determine the 
powers Xj t i for j = 1, . . . , k, and i = 1. . . . ,p. So the homogeneous system of linear equations has, in accord 



with Eq. (2.11 1, for each of these tt- invariants the alternative notation 



G • Xi = 



5i,i 

52,1 



01,2 
52,2 



5r,l 5r,2 



9l,K 

52, n 
9r,n 







II 









= {0} for£=l,.. 



(2.12) 



The rank of the dimensional matrix is equal to the number of fundamental dimensions, r. If the number of 
dimensional quantities, n, is equal to r, we will obtain: p = 0. In this case there is only a trivial solution. In 



the case of p > 0, the homogeneous system of linear equation (2.12) is indeterminate, i.e., more unknowns 
than equations, a fact that is true in all instances presented here. Hence, for each of the p non-dimensional 
7r numbers, it is necessary to make a reasonable choice for p of these unknowns, x Ki i, to put this set of 
equations into a solvable state. After that we obtain for each ir number an inhomogeneous linear equation 
system that serves to determine the remaining r — k — p unknowns. Thus, the remaining r x r dimensional 
matrix Go = {g r ,r} has the rank r, too. It is the largest square sub-matrix for which the determinant is 
unequal to zero (|g r .r| 7^ 0). Thus, we have 



Gn • X,; — Bj 



or 



51,1 


51,2 • 


• 5l,r 




Xl,i 


52,1 


52,2 • 


• 52,r 


II 


X2,i 


5r,l 


5r,2 • 


9r,r 







Bi, t 
B 2 ,i 



for i=l, 



(2.13) 



This inhomogeneous system of linear equations can be solved for x m ^ for m = 1, . . . , r by employing Cramer's 
rule. 



3 Similarity hypotheses for blackbody radiation laws 

In 1860 Kirchhoff 22J proposed his famous theorem that for any body at a given temperature the ratio of 
cmissivity, e„, i.e., the intensity of the emitted radiation at a given frequency v, and absorptivity, a v , of the 
radiation of the same frequency is the same generally expressed by 



(his 



(3.1) 
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where he called a body perfectly black when a v — 1 so that J(V, T) is the emissive power of such a black 
body An example of a perfectly black body is the 'Hohlraumstrahlung' that describes the radiation 

in a cavity bounded by any emitting and absorbing substances of uniform temperature which are opaque. 
The state of the thermal radiation which takes place in this cavity is entirely independent of the nature 
and properties of these substances and only depends on the absolute temperature, T, and the frequency (or 
the wavelength A or the angular frequency u> = 2 ir v). For this special case the radiation is homogeneous, 
isotropic and unpolarized so that we have |34j 

J(u,T) = U(u,T) . (3.2) 

Here, U(v,T) is called the monochromatic (or spectral) energy density of the radiation in the cavity. Kirch- 
hoff 's theorem has become one of the most general of radiation theory and expresses the existence of tem- 
perature equilibrium for radiation, as already pointed out by Wilhelm Wien during his Nobel Lecture given 
in 1911. Expressions for the monochromatic energy density called the blackbody radiation laws were derived 
by (a) Rayleigh [29j [30] and Jeans [20], (b) Wien [48] and Paschen [36], (c) Thiesen [46], (d) Rayleigh [29] . 
and (e) Planck [37] |M] 133- 



3.1 The Rayleigh- Jeans law 

The dimensional 7r-invariants analysis described in chapter 2 can be employed to derive all these blackbody 
radiation laws. First, we consider the radiation law of Rayleigh [29, 30 and Jeans 20J for which the similarity 
hypothesis F(Q 1 , Q 2 , Q3, Q4, Q5) = F{U,v,T,c,k) is postulated. Here, Qi — U, Q2 = v, Q3 — T, Qi = c, 
the velocity of light in vacuum, and Q5 = fc, the Boltzmann constant. Obviously, the number of dimensional 
quantities is n — 5. The dimensional matrix is given by 



-1 








1 


2 








1 





-1 


-1 


-1 





-1 


-2 


1 











1 



(3.3) 



that can be derived from the table of fundamental dimensions, 





U 


V 


T 


c 


k 


Length 


- 1 








1 


2 


Temperature 








1 





- 1 


Time 


- 1 


- 1 





- 1 


- 2 


Mass 


1 











1 



In accord with Eq. (2.12 1, the homogeneous system of linear equations is given by 



}={0} 



-1 








1 














1 









X2,l 


-1 


-1 





-1 




II 


#3,1 1 


1 















£4,1 
, X 5,l , 



(3.4) 



Obviously, the rank of the dimensional matrix is r — 4, and we have p = k — r — 1 non-dimensional 7r number 
that can be deduced from 



- Xx : i + x^i + 2 x 5i i = 

£3,1 - £5,1 = 

- X\s - X2,l - 2:4,1 - 2 £5,1 = 
£1,1 + £5,1 = 



(3.5) 
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This equation set results in X2,i — —2 11,1, 2:3,1 = —£1,1, X4.1 — 3 £1,1, and X5.1 = —£1,1. Choosing x\^\ = lQ 
yields (see Eq. Q) 



Rearranging provides finally 



7r 1 = nC' 1 = U l v- 2 ^ 1 c^k- 1 . (3.6) 
3=1 



U(v,T) = TT l ^rkT , (3.7) 



where the 7r-invariant can be identified as 7Ti = 87r. Thus, we have 

V (y iT )= % -^^kT . (3.8) 

This radiation law was first derived by Rayleigh [29, 30 using principles of classical statistics, with a correc- 
tion by Jeans [2D]. Lorentz [!JT] derived it in a somewhat different way. Obviously, it fulfills both Kirchhoff 's 
[22] findings and the requirements of Wien's [17] conventional displacement law, U{v,T~) oc i/ 3 f(u,T). Equa- 



tion (3.8 1 is called the classical blackbody radiation law because, as already pointed out by Sommerfeld 
[43j . it is restricted to the two constants c and k which are well known in classical physics. In the case 
when the frequency is small and the temperature relatively high, this formula works well (indeed it can be 
considered as an asymptotic solution for v — > 0), as experimentally proofed by Lummcr and Pringsheim [32 
and especially by Rubens and Kurlbaum [U] [32] . But it is obvious that the law of Rayleigh [2DJ [3D] and 
Jeans [2D] cannot be correct because for v — ► 00 the monochromatic energy density, U(v,T), would tend to 
infinity. Ehrcnfest [14] coined it the Rayleigh- Jeans catastrophe in the ultraviolet. Consequently, the integral 
over U(v, T) for all frequencies yielding the energy density, 

f°° 8 7T f°° 

E(T)= U{u,T) dv = — s-fe T I v 2 dv = oo , (3.9) 



would become divergent [15] , in complete contrast to Boltzmann's [Jj thermodynamic derivation of the T 4 
law, 

E(T) = aT 4 (3.10) 

and Stefan's [45] empirical finding 

R(T) = ^E(T) = ~ T 4 = a T 4 (3.11) 
where a = c a/4 is customarily called the Stefan constant. 

3.2 A generalized form of Wien's displacement law 

Since this classical blackbody radiation law completely fails in the case of high frequencies, it is necessary to 
look for an improved similarity hypothesis that must lead to more than one 7r-invariant so that a universal 
function (or similarity function) is established. Such a universal function is indispensable to ensure that the 



integral in Eq. ( 3.9 ) stays finite. However, we cannot introduce another independent variable to get more than 
one 7r-invariant because it would be in disagreement with Kirchhoff 's [22j findings. Since, on the other hand, 
the Raylcigh-Jeans law seems to be reasonable for small frequencies, one may state that the two dimensional 



) If we choose xi,\ = a, where q ^ is an arbitrary real number, we will lead to another invariant 7Tj. The relationship 



between this 7r-invariant and that, occurring in Eq. (3.6 I, is given by n\ = Therefore, for convenience, we may simply 

choose: Xi 1 = 1. 
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constants of the classical physics, c and k, should also occur in an improved similarity hypothesis. Conse- 
quently, we must look for another dimensional constant |43j and have to make a choice for its fundamental di- 
mensions. Thus, an improved similarity hypothesis reads F(Qi, Qi-, Q3, Qi, Q5, Qe) = F(U, v, T, c, k, rj) = 0, 
i.e., instead of the restriction to the two classical constants c and k, we now have three dimensional constants 
in our similarity hypothesis, where the dimensions of the third constant, Qq — rj, is to be chosen in such 
a sense that rj vj[k T) becomes non-dimensional, in hoping to establish another 7r-invariant, and, hence, a 
universal function ifR^rj vj(k J 1 )), according to the explicit representation tti — ^Pr{^2) expressed by Eq. 



(2.5). This means that not only k T, but also 77 v have the dimensions of energy. Obviously, the number of 
dimensional quantities is k, = 6, and the table of fundamental dimensions reads 





U 


V 


T 


c 


k 


X] 


Length 


- 1 








1 


2 


2 


Temperature 








1 





- 1 





Time 


- 1 


- 1 





- 1 


- 2 


- 1 


Mass 


1 











1 


1 



leading to the dimensional matrix: 



G = 




1 
-1 




(3.12) 



The homogeneous system of linear equations can then be expressed by (see Eq. ( 2.12| |) 



1 

1 

-1 -1 






>={0} fori = 1,2 



(3.13) 



where the rank of the dimensional matrix is r = 4, i.e., we have p — k — r = 2 non-dimensional 7r-invariants 
so that a universal function, as urgently required, is established. These two tt numbers can be derived from 



Xi 



X3,- 



for 



1,2 



(3.14) 



+ £4,; + 2 x b i + 2 xej 

^1.! X2 5 i X4 _i 2 £5.2 ^6 5 2 

Choosing x± t i — 1, X6,i = N, — 0, and 2^2 = 1 yields X2,i = —2 + N, 2:3^ = — 1 — N, x^ t \ = 3, 
x$ t i = — 1 — N, X2,2 — 1, X3 : 2 = — 1, £4.2 = 0, and £5.2 = — 1. Note that the choice x± t i = 1 and x\^ = 
is indispensable to ensure that U(y,T) only occurs explicitly. Furthermore, ./V < 3 is a real number, and 




x 6,2 — 1 is chosen in such a sense that 77 vj(k T) is established as a 7r-invariant. According to Eq. (2.2 1, the 
7r-invariants are then given by 



7Tl 



3=1 



and 



U 1 „-2 + N T -l-N J k -l-N 



U c 3 / ry v \ N 



v 2 kTXkT 



(3.15) 



n 2 = l[Q x /' 2 = U° v 1 T- 



V = 



•q v 

It 



(3.16) 
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According to Eq. (2.5 1, we have 



7Tl 



U c 3 



v 2 kT 

where <Pr(ji v /{k T)) is the universal function. Rearranging yields then 



U(v,T) 



v 



k T 



(r]v\- N fr)v\ 
\kT) * R \kT) 



For historical reasons and convenience, we may introduce 8 7r into Eq. (3.181. In doing so, we have 

r 2 

<■■• 



where 



ft] v 
l \kT 



8 7T 



\kTJ 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



The merit of Eq. (3.181 is that the monochromatic energy density U{y,T) is led to the universal function 
(pn(r) vj(k T)) which only depends on the non-dimensional argument 7r 2 = r\ vj(k T). Obviously, this 



equation substantially agrees with Kirchhoff's [25] findings. Equation (3.191 may also be considered as a 
generalized form of Wien's [17] displacement law. The conventional form of this law can be derived by setting 
N = — 1. In doing so, one obtains 



U(v,T) 



8 7T J] V 3 



4r) 



(3.21) 



3.2.1 Criteria for determining the universal function heuristically 

As mentioned before, the universal function ^r(t] vj(k T)) cannot explicitly be inferred from dimensional 
arguments. However, any form of such a universal function must be compatible with following criteria: (a) 
There must exist a distinguish maximum. (b)It must guarantee that the integral in 



E(T) = 8 7T k 



k 

c rj 



rp4 



X 2 - N <5> R (X) dX 



(3.22) 



with X = 7r 2 
identified as 



rj vj(k T) keeps finite. According to this equation, the constant a in Eq. (3.10) can be 



\crjJ 



X Z ~ N <$> R (X) dX 



(3.23) 



Note that Wien [48] also recognized the power law of Stefan [45] and Boltzmann [7J in deriving the quant ity 
ci v 3 that occurs in his radiation law. First, it must fulfill the requirement that the integral in Eq. (3.231 is 
convergent. If we assume, for instance, &r(X) = X N ~ 3 , then we will obtain for this integral 



/ X 2 - N <S> R (X) dX = lim / 

JO b ^ a SJa 



= lim In — 



X 



' oo 
- 



(3.24) 



Obviously, this integral is divergent, too. Consequently, the tendency of the integrand X 2 ~ n Qr(X) to 
zero when X approaches to infinity must be faster than X^ 1 . Ehrenfest [14], who only considered the case 
N = — 1, postulated that, therefore, the requirement 



lim A 4 $ R (X) = (3.25a) 
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must be fulfilled. He coined it the violet requirement (also called the violet condition |24j). For any real 
number of N < 3 we may slightly modify it to 



lim X^^^ji^^O . (3.25b) 

V — > oo 

A more rigorous requirement, of course, would be 

lim X m $ R (X) = , (3.26) 

V — ' oo 

where m > 3 — N. This means that for X — ► oo the universal function Q R (X) must tend with a higher 
intense to zero than X~ m for any m > 3 — N. For N = — 1, Ehrenfest [H] coined it the strengthened 
violet requirement, (c) The energy density, U{v,T), must tend to unity when v becomes smaller and smaller 
because the Rayleigh- Jeans law (3.8 1 is the asymptotic solution for that case. Following Ehrenfest [TJ], this 
requirement may be called the red requirement (also called the red condition |24j ) being expressed by 

.■aG?)"" •«(*?)-' • < 327 > 

3.2.2 The Maximum condition 



The first derivative of Eq. (3.19) reads 

U\v, T) = ^ ^~ N {(2 - N) <S> R (/3 v) + p v $ H '(/J u)} , (3.28) 

with a = c\ = 8 ir rj/c 3 and (3 = rj/(k T). Thus we obtain for the extreme of U(u, T) 

U\u e ,T)=0 «-> (2-N)$ R (pv e )+pv e $ R '(J3u e )=0 , (3.29) 
where v e is the frequency of the extreme. Defining x = /3 v e results in 

(2 - N) <f> R (x) + x §r\x) = . (3.30) 
The solution of the equation is given by 

<S> R {x)=x N - 2 . (3.31) 
For v e and a given temperature <& R (x) provides a certain value Z. Thus, one obtains 

v k i 

= - Z n-2 = const. (3.32) 
T rj 

that may reflect Wien's displacement relation z^ max /T = const, or A max T = const. Here, ^ max is the fre- 
quency, for which U(v,T) reaches its maximum, and A max is the corresponding wavelength. The second 
derivative reads 

IT (i/, T) = ^v- N {(1 -N)(2- N)<S> R (f] v) + 2 (2 - N) )3 v v) + {fi vf ((3 v)} . (3.33) 

Considering Eq. ( 3.31[ ) yields then for the extreme 

U"(y e ,T)=0 , (3.34) 

i.e., we have a stationary point of inflection, rather than a maximum (U" (v e ,T) < 0) as requested. This is 
valid for any finite value of N, i.e., power laws as expressed by Eq. (3.311 are excluded due to the maxi- 
mum condition. Since the maximum condition is not generally fulfilled, both Ehrenfest 's [14 statement that 
Wien's displacement law does not impose any restriction a priori on the form of <$> R (X) and Sommerfeld's 
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[43] derivation of Wicn's displacement relation are not entirely accurate. 

In contrast to Ehrenfest's statement, the condition U"(j/ e ,r) < clearly imposes a restriction on the 
form of Qr(X). For N < 2, for instance, the exponential function, 



Q R {X) =exp(-X) , 

guarantees that the maximum condition U"(y e ,T) < is always fulfilled (see Appendix B). 



(3.35) 



3.2.3 Ehrenfest's red and violet requirements and the blackbody radiation laws of Rayleigh, 
Wien and Paschen, Thiesen as well as Planck 



Using this exponential function (3.351, the integral in Eq. (3.23) becomes 



X 2 ~ N exp (-X) dX = T(3 - N) 



Here, Euler's T-function defined by 



r(x) 



P CO 

/ X x ~ l exp(-X) dX 
Jo 



(3.36) 



(3.37) 



for all real numbers \ — 3 — N > has been applied. This condition is clearly fulfilled when the restriction 
of the maximum condition, N < 2, is considered. Obviously, using this exponential function guarantees that 



the integral in Eq. (3.231 keeps finite. It also obeys the strengthened violet requirement. If we accept this 



exponential function for a moment, we will obtain 



U{v,T) 



8ttu 2 



k T 



\kr) exp \- It) 



(3.38) 



For N = 0, for instance, we obtain Rayleigh 's [29 radiation formula. In this case the integral in Eq. (3.361 
amounts to T(3) = 2. Choosing N = — 1 yields the radiation law of Wien 48J and Paschen 36J with 
r(4) = 6. In the case of Thiesen's (46] radiation law, which can be derived by setting N = — 1/2, we will 



obtain T(3.5) = 3.3234. Obviously, Eq. (3.381 contains the radiation laws of (a) Rayleigh, (b) Wien and 
Paschen, and (c) Thiesen as special cases. 

Apparently, both the Wien-Paschen radiation law and that of Thiesen obey the strengthened violet 
requirement and, of course, fulfill the maximum condition. However, they do not tend to the classical black- 
body radiation law of Rayleigh and Jeans given by Eq. (3.8 1, i.e., they do not obey the red requirement. 



If we namely express the exponential function by a Maclaurin series, we will obtain for the Wien-Paschen 
radiation law (N = — 1) 



X 2 X 3 

For small values of X or r\ v <C k T this series can be approximated by 



exp(X) = 1 + X 



(3.39) 



Thus, Eq. (3.381 may read 



exp(X) = 1 + X 



U{v,T) = -— k T 



1 + X 



(3.40) 



(3.41) 



Obviously, the expression X/{\ + X) does not converge to unity, as requested by Eq. (3.8) (see also Figure 
[J). The same is true in the case of Thiesen's radiation law that can similarly be approximated for that range 

by 



U(v,T) 



7T V 2 X^ 

kT 



1 + X 



(3.42) 
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v/Tin 10 11 s" 1 K" 1 



Figure 1: Planck's [331351135] function of the spectral energy density. The curve represents U(v, T)/T 3 versus 
v/T so that it becomes independent of T. Note that the constant, C — 5.8787 • 10 10 s _1 K -1 , occurring 
in Wicn's displacement relationship C = v msx jT can be inferred from this figure (dotted line). Also shown 
are the functions of Wien [JS], Paschen [35], Raylcigh [29 , as well as Rayleigh [25] |3U] and Jeans [2U] . The 
constants used for plotting these functions are listed in Table [T] 



On the contrary, if we choose 



and, again, N — — 1, Eq. (3.191 will provide 



exp (A) - 1 



(3.43) 



U{y,T) 



7T rj v 



exp 



rj v 
If 



1 



(3.44) 



This equation is quite similar to Planck's [371 EFJ [39] radiation law. If we again express the exponential 
function by a Maclaurin series and again assume rj v <^ k T, we will obtain 



U(v,T) 



k T 



X 



8ttv 2 



1+X-l 



k T 



(3.45) 



Apparently, Eq. (3.441 completely fulfils the red requirement (sec Figure]!]). For rj v kT ', i.e., exp (X) 1, 
one also obtains the Wien-Paschen radiation law. This fact simply states that Eq. ( |3.44[ ) also fulfils the 
strengthened violet requirement (see Figure [T]). 



Now, we have to check whether the integral in Eq. (3.23 ) keeps finite. Applying &r(X) = (exp (X) — 1) 
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Table 1 : Quantities used for plotting the functions illustrated in Figures [T] and [2] 



Author 


k 

(J if" 1 ) 


V 

(J s) 


C 

(s- 1 K- 1 ) 


Planck 


1.3806 ■ 10- 23 


6.6262 ■ 10- 34 


5.8787 • 10 10 


Wien — Paschen 


1.7963 • 10- 23 


9.1670 • 10- 34 


same 


Thiesen 


1.8768 • 10- 23 


7.9813 • 10- 34 


same 


Rayleigh 


1.5967- 10- 23 


5.4323 • 10- 34 


same 


Rayleigh — Jeans 


1.3806 • 10- 23 


6.6262 • 10- 34 





and assuming N — — 1 yield then 



V dX = r^exp(-X) x 



(3.46) 



exp (X) - 1 J 1 - exp (-X) 

poo 

= / X 3 exp(-A)(l+exp(-A) + exp(-2A) + ...)dA 
Jo 

= J X 3 (exp(-X) + exp(-2X)+exp(-3X) + ...) dx} 
Substituting nX (n = 1, 2, 3, . . .) by Y yields then 

where the sum in this equation 

00 1 2 2 k - 1 7T 2 fc 



can be calculated using the Bernoulli number for k = 2, namely B 2 = 1/30. Apparently, the integral (3.471 
is convergent. As a and, hence, a = 4 a/c are known, the constant n can be determined. In doing so, 
one obtains: n — 6.6262 • lO -34 J s, i.e., it equals the Planck constant h. Note that Planck called it the 
"Wirkungsquantum" that means an elementary quantum of action, and the product h v is customarily des- 
ignated as the quantum of energy. Furthermore, the expression 

(3.49) 



that occurs in Eq. (3.441 is customarily called the Planck distribution. It may be considered as a special case 
of the Bose-Einstein-distribution when the chemical potential of a "gas" of photons is considered as u = 

nana 1231 ■ 
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10.0 



Planck 

Thiesen 

Rayleigh-Jeans 




v/Tin 10 11 s" 1 K" 1 



Figure 2: As in Figure [T] but the functions of Wien 
Thiesen \M- 



Paschen [33], Rayleigh [25] are replaced by that of 



It is obvious that Planck's [37J EB 39J law obeys the requirements in the red range and in the violet range, 
and, in addition, it fulfils the maximum conditions. The same, of course, is also true in the case of Rayleigh's 
[29) radiation formula; but as illustrated in Figure [I] there are appreciable differences between both radiation 
laws in the higher-frequency range. These differences can also be inferred to the value of the integral in Eq. 
(3.361; for N — 0, it amounts to T(3) = 2. In that range also notable differences exist between Planck's 
law and those of Wien [JS] and Paschen [SB] because of similar reasons (see Figure [l}. Since 7r 4 /15 = 6.49, 
Planck's distribution yields a value of the integral in Eq. (3.231 which is notably larger than that provided 



by this equation when for the asymptotic solution for h v 3> k T ', namely the radiation law of Wien and 
Paschen, N = — 1, is chosen. As illustrated in Figure [2] the differences between the radiation laws of Planck 
and Thiesen [IB] are appreciably smaller than in the other cases. These differences are so small that empirical 
results alone may be insufficient for evaluating both radiation laws. Since Thiesen's radiation law, of course, 
does not fulfill the red requirement, Planck's one has to be preferred. Even though the radiation law of Wien 
and Paschen represents an asymptotic solution for Planck's law, the former hardly represents a better solution 
than Thiesen's one. One may speculate that the small differences between the radiation laws of Planck and 
Thiesen would be ignored today. Note that the functions of Planck and Rayleigh-Jeans illustrated in Figure 
[TJare based on the values of k = 1.3806 • 10~ 23 J K^ 1 and r\ = h = 6.6262 • 10~ 34 J s currently recommended. 
The corresponding values for the functions of Wien, Thiesen, and Rayleigh have been derived using Stefan's 
constant a = 5.6696 • 1CT 8 JmT 2 s^ 1 R- i and Wien's relation v max /T = const. = 5.8787 • 10 10 s" 1 K~ x . 
This constant was calculated using X = 2.82144, iteratively computed on the basis of the transcendental 
equation X = 3 (1 — exp(— X)) obtained from Planck's radiation law as the condition for which U(y,T) 
reaches its maximum 128 , and the recommend values of k and h. 
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4 Some historical notes 

Planck presented his blackbody radiation law at a meeting of the German Physical Society on October 19, 
1900 in the form of [37J HO dD 131 HO] , 

U(u,T) = Cl % (4.1) 

exp(c 2 -) - 1 

On the contrary Eq. (3.44) reflects the form of Planck's radiation lawj^Jas presented in his seminal paper 
published at the beginning of 1901 [3H] and presented to the German Physical Society on December 14, 1900 
[33 O [27l I3H HQ]. This day may be designated the birthday of quantum theory because the elementary 
quantum of action explicitly occurred in Planck's radiation law (e.g., |27[ 134] . Our use of principles of 



dimensional analysis in heuristically deriving Eq. (3.441 by ignoring the aid of the linear harmonic oscillator 



model and Planck's assumption that the energy occurring in Boltzmann's [5] distribution is quantized gives 
evidence that Planck's findings were the results of a lucky chance. In his Nobel Lecture, delivered in 1920, 
Planck objectively stated: ". . . even if the radiation formula should prove itself to be absolutely accurate, 
it would still only have, within the significance of a happily chosen interpolation formula, a strictly limited 
value." 

The first who, indeed, realized the true nature of Planck's constant was Einstein [TS]. In his article he 
related the monochromatic radiation, from a thermodynamic point of view, to mutually independent light 
quanta (or photons) and their magnitude, in principle, to Ae = h v that occurs in Planck's radiation law 



(see Eq. ( 3.44 1) . Einstein, however, did not deal with Planck's radiation law in its exact manner, but 
rather with the approximation that fits the radiation law of Wien [48] and Paschen [36]. It seems that 
he first recognized that the quantum discontinuity was an essential part of Planck's radiation theory (e.g., 
[TTl [231 |2"T] ) . As discussed by Klein [22] and Navarro and Perez [33] , a milestone on the road of quantum 
discontinuity and light quanta is Ehrenfest's [H] article on the essential nature of the different quantum 
hypotheses in radiation theory. With his article Ehrenfest contributed to the clarification of the hypothesis 
of light quanta. Unfortunately his contribution was not recognized for a long time. 



Appendix A: Jeans' attempt to derive Wien's displacement law using dimen- 
sional analysis 

Following Jeans [HI [21], the similarity hypothesis reads F(Qi, Q2, Q3, Q4, Q5, Qe, Q7, Qs) = F(U, X,T,c,e, 
m,R,K) = 0. Here, Q2 = A is the wavelength, Q 5 = e and Qq = m are the charge and the mass of an 
electron, respectively, = R is the universal gas constant, and Qs = K the dielectric constant of the 
ether expressed with respect to an arbitrary measuring system. All other symbols have the same meaning as 
mentioned before. Since the monochromatic energy density has to be considered, the number of dimensional 
quantities is k = 8 now. The dimensional matrix is given by 



G 



3/2 


-1 
1/2 
1/2 



-1 
-2 

1 





(Al) 



that can be inferred from the table of fundamental dimensions given by 



2 ) Note that Planck derived Eq. j3.44| using an equation for the entropy S of a linear harmonic oscillator, and he related it 
to its mean energy and to the quantum of energy, h v. 
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U 


A 


T 


c 


e 


m 


R 


K 


Length 


- 1 


1 





1 


3/2 





2 





Temperature 








1 











- 1 





Time 


- 1 








- 1 


- 1 





- 2 





Mass 


1 











1/2 


1 


1 





Arbitrary system 














1/2 








1 



The homogeneous system of linear equations can then be written as (see Eq. (2.12 1) 



-1 1 1 3/2 2 ' 

10 0-10 
-1 -1 -1 -2 I < 

1 1/2 1 1 
1/2 1 



Xl,i 
%2,i 

X^,i 

XQA 

X S .i 



{0} fori = 1,2, 3 



(A2) 



Since the rank of the dimensional matrix is r = 5, we have p = k — r = 3 non-dimensional 7r-invariants. 
Thus, a universal function of the form 7Ti = <p(jr 2 , "3) is established. These three 7r-invariants can be derived 
from 



- xi t i + x 2 ,i + x^ + 3/2 x 5>i + 2 x 7 ,i 

%3,i - X7,i 

X%,i X^ i X§^ 2 X*j^ 

X\,i + 1/2 x 5ii + X 6 ,i + X 7 ,i 

1/2 x 5li 



X8, 



= 

= 

= 

= 

= 



for £ = 1,2,3 



(A3) 



Choosing x M = 1, x 6)1 = 0, x 7>1 = -1, x 1(2 = 0, x 6a = -1, £7,2 = 1, ^3,1 = 0, ^6,3 = ~ 2 and 2:7,3 = 1 
yields x 2 ,i = 2, £3.1 = -1, 2:4,1 = 1, £5,1 = 0, x a ,i = 0, 2:2,2 = 0, 0:3,2 = 1, 2:4,2 = -2, 2:5,2 = 0, x s , 2 = 0, 
2:2,3 = —1, 2:3.3 = 1, 2:4,3 = —4, 2:5 3 = 2, and 2:8,3 = — 1. The choice 2:1,1 = 1, 2:1.2 = 0, and x\$ — is re- 
quired to guarantee that U (A, T) only occurs explicitly. In accord with Eq. ( 2.2 1 , the 7r-invariants are given by 



m = J] Q\ = U 1 A 2 T~ 1 c 1 e° m° R~ 1 K° = 

8 

7T 2 = JJ Q] " = U° A° T 1 c- 2 e° m- 1 R 1 K° 



l/A 2 c 
RT 

RT 



m cr 



(A4) 
(A5) 



and 



fiTe 2 
A m 2 c 4 if 



According to Eq. (2.5 1, we have 



7Tl 



t/A 2 c 
RT 



RT RT e 2 



m c 2 ' A m 2 c 4 if . 



(A6) 



(A7) 



or 



(A8) 
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It is apparent that Eq. ( A8 1 completely disagrees with 



f/(A,T)cx~^(A,T) (A9) 



U(\, T) cx A~ 5 ipr{X, T) (A10) 



that can be derived from Eqs. (3.191 and ( |3.21 I, respectively. Since we always have Xi^ = 2 Xxj, any other 
choice can give no better relationships. Therefore, Jeans' attempt to derive Wien's jlTj displacement law 
completely fails because of his inadequate similarity hypothesis. 

Appendix B: Derivation of Wien's displacement relation for = exp (—X) 

If we consider the exponential function v) = exp (— [3 v), the first two derivatives of Eq. ( |3.19[ ) will read 

U'(u, T) = j„ (2-N-fiv) exp (-/3 v) (Bl) 

and 

U»( Vl T) = ^v- N {(\-N)(2-N)-2(2-N)Pv + ((3v) 2 } exp (-/? z/) . (B2) 
For the extreme we obtain [3 v e — 2 ~ N ov 

^ = (2-N)- = const. (B3) 
T rj 



Introducing Wien's displacement relation into Eq. (B2| yields finally 



IT {v e , T) = (N - 2) exp (N ~ 2) , (B4) 

This means that U"(iy e ,T) becomes negative for any real value of N that fulfils the condition N < 2, i.e., in 
such case the condition of a maximum is fulfilled. 
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